Effects of paramagnetic impurities on two-band superconductor 
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I calculated the effect of magnetic impurities on the normal and superconductive properties of a 
multiband s-wave superconductor by direct solution of the two-band Eliashberg equations. In this 
way I determined the critical temperature, the values of the superconductive gaps, the shape of the 
superconductive density of states and other physical quantities that depend on the concentration of 
magnetic impurities. I found that the gaps and the penetration lengths display an unusual behaviour 
as a function of temperature. I examine the possibility that the presence of a negative induced gap 
raises the critical temperature. 
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A specific, but representative case of anisotropic su- 
perconductivity is multiband superconductivity where 
the order parameter is different in different bands. The 
more famous and clear case is the magnesium diboride [l| 
where the two gaps are observed in several experiments. 
Some properties of this material changes markedly from 
the BCS s-wave one band case. For example the rate 
between the gaps and the critical temperature can be 
major or minor then 3.53. Another difference is that the 
non magnetic interband scattering rate produces a reduc- 
tion of the critical temperature [2| while in an one band 
superconductor there isn't effect on the critical temper- 
ature. Recently also the effects of magnetic impurities 
on a multiband superconductor have been studied, in a 
theoreticaljl, 0] and experimental way, on MgB2 doped 
by Mn 0,[|7 

Let us start from the generalization of the Eliashberg 
theory Q for systems with two bands 0, Q ,that has al- 
ready been used with success to study the MgB2 system 
S EH EH E2 • To obtain the gaps and the critical tem- 
perature within the s-wave, two-band Eliashberg model 
one has to solve four coupled integral equations for the 
gaps Ai(iuj n ) and the renormalization functions Zi(iu n ), 
where i is a band index and u>„ are the Matsubara fre- 
quencies. I included in the equations the non-magnetic 
and magnetic impurity scattering rates in the Born ap- 
proximation, and Tfj. 

u> n Zi (iuj n ) = 0J n + ttT 2J Ajj (iuj n - iu) m )N j z (iu) m ) + 

m,j 
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Zi(iu! n )Ai(iuj n ) = itTyi[Aij(iu n - iu m ) - t4i( u c)] • 
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where 9 is the Heaviside function, oj c is a cut-off energy 
and Aij(iu n - iu m ) = f °° c^a|--F(w)/[(a;„ - uj m ) 2 + 



lu 2 ], N 3 A (iu m ) = Aj(iw m )/^/w^ + A 2 (iuj m ), N z (icu m ) = 

w m / yjuj^ + A 2 (iuj m )- 

The solution of the Eliashberg equations requires as in- 
put: i) the four (but only three independent pj) electron- 
phonon spectral functions a?-(w)i f, (w); ii) the four (but 
only three independent 0]) elements of the Coulomb 
pseudopotential matrix /i*(w c ); iii) the two (but only 
one effective Q) non- magnetic impurity scattering rates 
T^; iv) the four (but only three independent Q) para- 
magnetic impurity scattering rates T^. The four spec- 
tral functions a 2 Ju>)F(u>), that were calculated for pure 
MgBz in ref. [ll|, have the following electron-phonon 
coupling constant: A crcr (x = 0)=1.017, A 7rTr (a; = 0)=0.448, 
A CT7r (x=0)=0.213 and X w<T (x=Q)=0.156. 

As far as the Coulomb pseudopotential is concerned, I 
use the expression calculated for pure MgB2 fl3j |: 
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where n{ui c ,x) is a free parameter and N^ >t (E F ,x) is 
the total normal density of states at the Fermi level. For 
obtaining the experimental critical temperature of pure 
MgB 2 case (T c = 39.4 K) we fix fi(u c ) = 0.03105 with 
cut-off energy lu c = 500 meV and maximum energy 570 
meV. In all our calculations I use [TO] N$(E F ) = 0.30 
states/ (cell eV) and Nfi(E F ) = 0.41 states/ (cell eV). In 
this work I study only the effect of paramagnetic impuri- 
ties so is always = 0. In fig. 1 we can see the effect of 
magnetic impurities on the critical temperature in four 
limit cases. The stronger reduction of T c is when the 
is preponderant while the effect of the is week. 
A part the cases a small amount of impurities is safe 
for reduce T c in a considerable way. This is consistent 
also with the use of the Born approximation and with 
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FIG. 1: The calculated critical temperature as function of T^_: 
/ (solid line), TV / (dashed line), Y1 a = / 
(dotted line) and T1 n (dashed-dotted line). 

the fact of thinking the MgB% doped with paramagnetic 
impurities as a system perturbed and not as new ma- 
terial with different electron-phonon coupling constant, 
phonon energies and Coulomb pseudopotentials. In fig. 
2 we can see the calculated penetration depth by solution 
of imaginary axis Eliashberg equations [14j . The curve 
are remarkably different of pure case especially when is 
present interband scattering. 

In fig. 3 panel a we can see the calculated values of 
Ai(iuj n= o) as function of the temperature always with 
T c = 28K. When are present magnetic impurities in 
the interband channel the curve has a maximum more 
evident if the impurity content is bigger. In the ^ 
case the A w (ioj n= o) gap is negative in a small range of 
temperature (experimentally it is found that the presence 
of magnetic interband scattering is very small [1]). In the 
panel b we can see the calculated tunneling conductance 
at T = 4K: in this case the curves are similar for the 
effect of smearing of the temperature: it is necessary to go 
to lower temperatures for noting appreciable differences. 
Of course from experiments of quasiparticle tunneling it 
is impossible to determine the sign of A^. In the fig. 
4 it is possible see the calculated values of Ai(iui n= o) as 
function of the temperature with T c = 18K (panel a) and 
T c = 13.3K (panel b). In the panel b we can note that, 
in the ^ case, there is a range of temperatures 
where A 7r (iw„ =0 ) > A cr (iw„ =0 ). 

It is important remember that in presence of strong 
coupling electron-phonon interaction or magnetic impuri- 
ties the value of Ai(iu n= o) obtained by solving the imag- 
inary axis Eliashberg equations and the A^ obtained by 
real axis formulation can be very different so it is nec- 
essary always or the analitical continuation [l5j or the 
real axis solution for determine the physical quantities 
comparable to measurable experimental gaps. 

If we put the interband electron phonon coupling con- 
stant Acttt (of course also A^o-) arid Coulomb pseudopo- 
tential matrix fj,*^ equal to zero, in the MgBi case, I find 
an hypothetic material with T c = 45.55 K and only a 
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FIG. 2: The calculated penetration depth as function of tem- 
perature: panel a: pure MgB^ case (solid line ab plane, 
dashed line c axis), TZf 7^ (dotted line ab plane, dotted- 
dashed c axis); panel b: pure MgBi case (solid line ab plane, 
dashed line c axis), W = W 7^ (dotted line ab plane, 
dotted-dashed c axis);panel c: pure MgBi case (solid line ab 
plane, dashed line c axis), TV 7^ (dotted line ab plane, 
dotted-dashed c axis). 



superconductive gap in the a band because the 7r band is 
superconductive only for the presence of interband term. 
In the MgB2 case, the presence of interband term reduces 
T c . If I calculate the critical temperature in function of 
Xij with all other parameters equal to MgB% case, I found 
that when A^ = Xf$ B ' 2 the T c is next to a minimum as 
we can see in fig. 5. 

Now, if is only Ajj = 0, i 7^ j but the matrix fi*j 
has the usual MgBi values it can induce a negative gap 
in the 7r band and overall raise the critical tempera- 
ture until 48.75K. We can think of obtaining a negative 
A w (?o-' I i = o) gap i.e., roughly, when applies the condition 
Xij — [i\j < 0, i 7^ j, with chemical or field effect dop- 
ing where x is the doping content. I assume, for sim- 
plicity, X nu (x) = \ wa (0)N$(E F ,x)/N$(E F ,0) and 
AC(x) = /C(0) • N™(E F ,x) ■ NZ(E F ,0)/(N£(E Fl x) ■ 
N^ ot (E F ,0)) so I find, in the MgB 2 case, N£{E F ,x) < 
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FIG. 3: Panel a: the calculated Aj(ioVi=o), when T c — 28K, 
as function of temperature: YZ" (solid line), YZ a = / 
(dashed line) , YZJ 7^ (dotted line) ; panel b: the calculated 
tunneling conductance, at T — 4K, in the three different 
magnetic doping case: YZ" (solid line),r°: CT = TZT 
(dashed line), YZ? (dotted line) 
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FIG. 4: Panel a: the calculated Ai(iio n= o), when T c — 18K, 
as function of temperature: YZ" =£ (solid line), YZ° = YZ™ ± 
(dashed line), YZT 7^ (dotted line); panel b: the calculated 
Ai(iu n= o), when T c — 13.3K, as function of temperature: 
YZ" j= (solid line), YZ" = YZJ (dashed line), Y°7 ± 
(dotted line)). 



0.0311N$(E F ,x)/(0.3805N$(E F ,x) - 0.0311). For ex- 
ample in the case of Al doping, where, of course, I ne- 
glect the presence of disorder, i.e. I put = 0, this 
condition is safe only for x > 0.46 when the material 
isn't more a superconductor! To raise the critical tem- 
perature is more difficult because T c decreases princi- 
pally if \cc lowers and, in our rough model, \ aa {x) = 
\ aa .(0)N^(E F ,x)/N§(E F ,0) decreases with N$(E F ,x) 
so it is necessary to study in depth the problem. An 
other way for having a negative gap is if roughly \j — 
— Ffj < 0, i 7^ j, but the critical temperature, in this 
case, decreases: i.e. it is possible to have an induced neg- 
ative gap or with interband magnetic impurities or when 
the interband coupling is negative for example without 
phononic components but with the interband Coulomb 
pscuclopotcntial different from zero. The difference in 
the two cases is in the structure of Eliashberg equations: 
the Coulomb term is present only in the equation of the 
order parameter (eq. 2) and a negative gap in the channel 
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FIG. 5: The theoretical critical temperature of an hypothet- 
ical superconductor with all parameters less A CT7 r equals to 
MgB2 case. 



4 



7r produces a positive contribute in the more important 
channel a so the critical temperature raises while in the 
case of interband impurities in the equation of the order 
parameter there is a positive contribute of the Coulomb 
term but a negative contribute of the impurity term and 
also, in the equation of renormalization fuction (cq. 1), 
a positive contribute of the impurity term that raises 
the value of the renormalization function but, for conse- 
quence, lowers the critical temperature. 

I want remember that the possibility of a negative in- 



duced gap in a multiband system have been already pre- 
viewed by A. A. Golubov et al. more ten years ago (sec 
ref. 8). 

In conclusion the effects of magnetic impurities can 
produce, in multiband superconductors, unusual be- 
haviour as the temperature dependence of the order pa- 
rameter and penetration depth but overall can induce 
superconductivity with negative order parameter in the 
eventual normal band. 
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